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Abstract 

5 = 3/2 system with general isotropic nearest-neighbor exchange within a mean-field approxi- 
mation possesses a magnetically ordered ferromagnetic state and antiferromagnetic state, and two 
different spin nematic states, with zero spin expectation values. Both spin nematic phases display 
complicated symmetry break, including standard rotational break described by the vector-director 
u and specific symmetry break with respect to the time reversal. The break of time reversal is 
determined by non-trivial quantum averages cubic over the spin components and can be described 
by unit "pseudospin" vector a. The vectors a on different sites are parallel for a nematic state, 
and <r's are antiparallel for different sublattices for an antinematic phase. 

PACS numbers: 75.10.Jm, 75.40. Cx, 75.40.Gb 
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Quantum spin systems have provided a wide playground for the quest of novel types 
of quantum ordered states and quantum phase transitions. A number of exotic states has 
been discovered, the most well-known examples are the famous Haldane phase in integer 
spin antiferromagnetic chains [l| and quantum spin nematic states in the spin-1 bilinear- 
biquadratic isotropic magnets. A quantum spin nematic state has zero value of dipolar 
spin order parameter, (5) = 0, but spin rotation symmetry is spontaneously broken due to 
nontrivial quadrupolar spin expectation values of the type of (SiSj + SjSi), i,j = x, y, z. For 
5 = 1 spin operators, the quadrupolar variable have uniaxial symmetry and can be written 
through the condition ((5, w) 2 ) = 0. The order parameter can be present via the vector 
director u, in complete analogy to nematic liquid crystals. In the last two decades, this state 



has been actively studied for description of crysta 



low-dimensional magnets, see recent articles [4], |], [f| . The interest has got a considerable 
impact, motivated by investigation of multicomponent Bose-Einstein condensates of neutral 
atoms with nonzero integer spins 

aaa. 

The investigation of spin nematic states for 
5=1 systems has been performed for two different physical models, spin lattice system and 
and Bose gas of interacting particles with non-zero spin. The results obtained within both 
approaches are in a good agreement and complement each other. 

A few novel spin nematic states are found for Bose gas of interacting atoms with spin 



line magnets, see for review 

aa 

, including 



5 = 2, including the state with non- uniaxial symmetry [lOj. The question of the possibility 
of some kind of spin nematic states for half-integer spins (the minimal half-integer spin value 
allowing such state in 5 = 3/2) is of large general importance for physics of Fermi systems. 



As was found in Refs. 



11 



12 



13L high spin fermionic systems (ultracold Fermi gases) in some 



parameter region also exhibit many properties common to that for spin nematic states. On 
the best of our knowledge, for higher spins 5 > 1, both integer and half-integer, the only 
approach based on the direct analysis of gas of interacting particle has been used. 

The aim of this Letter is to develop a mean-field analysis of the ground state of spin 
5 = 3/2 isotropic spin Hamiltonian, with special attention to the spin nematic state, and 
to discuss the symmetry and the excitation spectra for such states. 

The most general isotropic exchange interaction for pare of spins 5 = 3/2 includes, 
additionally to Heisenberg bilinear interaction, non-heisenberg (biquadratic and bicubic) 
terms as well, which naturally leads to the model described by the following nearest-neighbor 
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interaction Hamiltonian: 
H = 



-E 

(0) 



+ L ( SfSfo 



(1) 



where Sj-are spin-3/2 operators at the lattice site I; J, K, and L are the exchange constants, 
corresponding to the bilinear, the biquadratic, and the bicubic exchange interactions, re- 
spectively, and summation over pares of the nearest neighbors 5 is implied. 

The spin-3/2 state \ipi) is a superposition of four basis states S z \tp s ) = s\ip s ), with 
coefficients z s , s = 3/2, 1/2, —1/2 and —3/2. The essential quantities are the ratios of 
these complex numbers, and the system parameter manifold is three-dimensional complex 
projective space CP 3 . It is convenient to present the ratios of z s as, 



-3/2 



ZZ/2 
Zl/2 
Z3/2 



vp 



tan 



e 

2' 



1/2 



Zl/2 



ip 



tan 



a 



2 ' 



iv 4 a 9 

3 ' • tan a ■ cos — • sec - 
^22 



(2) 



It is easy to check that the resolution of identity / D[ip] = 1, with the proper 

measure D[ip], is satisfied, and this six-parameter state has the properties of SU(4) coherent 
states. 

, one can construct the coherent state path integral and find the 



Using the states (|2J 
effective Lagrangian £ 



C = h J2\ \dtl + (d t p) sin 2 ^ 



+ (dt<p) sin 2 - cos 2 \i 



W, W 



sin /j, + 



H 



(3) 



where d t = d/dt and W is the "classical" (mean-field) energy of the system, which equals 
to the quantum mean value of the Hamiltonian with the states (j5J). 

The analysis of the ground state can be considerably simplified by making a rotation 
in the spin space to the principal axes of the spin-quadrupolar ellipsoid, i.e. by the use of 
condition (S l S j + S^S 1 ) — 0, if i ^ j. This condition yields the relations 9 + a = n, (3 = —7, 
ip = 7, and we can operate with the three-parameter state 



COS jl 



a 



sin 



+ e ' cos • 



a 



+ 



+ sin fj, 



. 01 
sin — 

2 



\ + e i7 cos^ 



(4) 



3 



2~ 



FIG. 1: (Color online) The domains of stability for different phases present through "angular 
variables" (fTTj) . 

For this state, diagonal biquadratic expectation values are independent on 7 and a, they 
can be written as 

4 (Si) = 3 + 2v / 3sin2/i + 4sinV, 

4<S 2 ) = 3 - 2V3 sin 2/^ + 4 sin V (5) 
The mean values of the spin operators are determined by all three parameters, 



(S x + iS y ) = sin n sin a- 



e 11 sin fj, + e n \/3 cos \x 



2(S Z ) = cos a [l — 4 cos 2 //] 



(6) 



and the angular variable 7 determines the orientation of the mean value of spin operator 
(S) in the xy— plane. 

Assuming a uniform ground state, one arrives to the following expression for the energy 
of the system W 

/ / • 1 no \ r ■> 
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(7) 



In fact, Eq. ((7|) gives us the free energy of the non-heisenberg S — 3/2 magnet with the 
Hamiltonian ([!]) at zero temperature T. Minimizing W over parameters /1, a and 7, one can 



find possible states of the system. Our analysis shows, at Ai > 0, where 



Ai = 16 J — 8K + 103L , (8) 

the minimum corresponds to the ferromagnetic phase with the saturated spin. For this state 
a — 7T, fM — 0, and 



| ^ferromagnet ) 



~), OT = |. (S x ) = (S y ) = 0. (9) 



Otherwise, at Ai < the minimum corresponds to the spin nematic state with a = tt/2, 
[A — 0, (S) = 0, and 



I ^ nematic ) \/2 



<(s 2 ) 2 ) = |, ((^) 2 ) = <(^) 2 > = |. (10) 

The image of this state is an elongated ellipsoid of evolution which can be described by 
unit vector-director u, u is parallel to z— axis for ffTUl) . So, the quantum phase transition 
"ferromagnetic state - spin nematic state" occurs at Ai = 0. In a natural parametrization 
of the exchange parameters through "angular variables" of the form 

J = JcosO, K = J sin 6 cos $, L = JsinOsin$, (11) 

the set of system parameters can be presented as a point on a sphere J 2 + K 2 + L 2 = J 2 . The 
condition Ai = is presented by a big circle on the sphere, passing through the direction 
p with tanO = 2J/K and $ = 0, and the point with L = 8K/103 on the equator, see 
Fig. 1. On this line, all uniform states of the system are continuously degenerated on the 
parameters a and /i. Thus there are lines with high symmetry in the parameter space of 
the S = 3/2 magnet. Note the equality Ai = corresponds to the condition of lifting of the 



. 3. 



symmetry of spin S = 3/2 Fermi gas in one band Habbard model, see Eq.(55) of Ref. 

The complete investigation of the stability of any phase with respect to arbitrary 
small perturbation should include an analysis of evolution of small deviations from the 
uniform state, Sa r = ot-r— o^o, where denotes one of six parameters of the state 
(J2J), oiifi corresponds to the ground state. Following the Bloch theorem, one can write 
8a.f— £g[5a:£ • exp(ifcZ)], k takes values within the first Brillouin zone. Information about 
stability can be extracted from the spectrum of elementary excitations to = Ui(k) found 
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for the state of interest. Such excitations can be learned by use of different approximate 



quantum approaches, e.g. bosonization of spin operators with the 1/N expansion 14[, or the 



Habbard operator technique, see, for example, 15|. Alternatively, one can apply a semiclas 



sical treatment based on the Lagrangian ([3]), with the usage of the concrete form of W, that 



is consistent with mean-field ground state calculations We had used two last methods 
for a d— dimensional hypercubic lattice, and found (as for spin S = 1 magnets, see [if]]) the 
same results for the both approaches. 

Variation of the Lagrangian ([3]) leads to the Hamilton system for generalized coordinates 
^7fc> an d its solution gives three branches of elementary excitations. Here we discuss 
briefly the properties of modes important for understanding of the stability of states. 

For the spin nematic state, the properties of two spin wave modes are common to that for 
a spin S — 1 nematic, their spectra are degenerated with the gapless dispersion law, linear 
over k = \k\ at k — > 0, uji^(k) — > ck, where he = (3a/8)^—2zAi ■ (AK — 5L). The third 
mode also has a gapless spectrum with linear asimptotics at k — > 0, 



= \^Jh HAi + 16 J £ - 8^ + 119^) , (12) 

where z is a coordination number, Jr, Kr and are Fourier-components of exchange 
constants of the form 

J k = J J2 S (} - e<ss ) ' J k - J ( ak ) 2 at k • 

The phase speed C3 of this mode vanishes at L — ► 0, Hc^ = (9a/4)y/—zLAi. For A x < 
these modes describe a long- wave instability at the line (jSJ). Additionally, at A 2 < 0, where 

A 2 = 16 J — 8if + 135L , (13) 

the values of 013 are imaginary for the quasimomentum k at the edge of Brillouin zone, where 
Jfc = 2zJ, ect. This means the instability, which leads to the transition to two-sublattice 
phase (antiferromagnetic, see below) on the line (fTBl . Thus, the spin nematic state is stable 
if and only under the conditions A x < and A 2 > 0. Note the line A 2 = describes a big 
circle, passing through the same direction p as for Ai = 0, see Fig. 1. Thus, the nematic 
state possesses three Goldstone modes (not two as for S = 1 system). We will discuss this 
fact below. 

The situation for the ferromagnetic state is simpler. One mode corresponds with the 
spin oscillations with \(S)\ = 3/2; it have the gapless dispersion law, parabolic at k — » 0, 



16hui(k) = (ak) 2 (48 J + 72K + 189 L) , and standard for Heisenberg ferromagnets. Two other 
modes possess oscillations of 1(5)1. They have gaps, proportional to Ai and describe the 
long-wave instability at Ai < 0. As well, these modes at A 2 < show the instability with 
respect to the transition to a two-sublattice state (an antinematic phase, see below). Thus, 
ferromagnetic state is stable if and only A x > and A 2 > 0. 

Let us return to spin nematic state. An appearance of the third Goldstone mode (fl2l) for 
this state can not be explained within a standard imaging of a nematic state as a ellipsoid 
of rotation, connected with bilinear spin expectation values (flOl) . But for spin 5 = 3/2 
magnet, an additional break of symmetry is caused by non-trivial expectation values cubic 
over spin components. For the state (jl]), such values are ov + ) = (1/3)((S X + iS y ) 3 } and 
<7(_) = (1/3) ((S x — iS y ) 3 }, <7(±) = cos 2 jj sin a exp (±ry). The transformation properties of 
the quantities a x = (c(+) + <7(_))/2 and a y = — cr(+))/2 are the same as for the 

components of the planar (two-dimensional) vector a, a oc e x cos 7 + e y sin 7 under the 
rotation around z— axis, \a\ = 1 for the nematic state. Having in mind that a changes its 
sign at time reversal t — > —t, we can say that the properties of this unit vector are common 
to that for spin expectation value and the mode ( 1T21) is associated with the oscillations of 
this vector. 

The presence of such characteristic as a is a feature of principal importance. As we 
mentioned above, for a spin S = 1 nematic, the symmetry breaking is associated with the 
quadrupolar ellipsoid only. But the spin 5 = 3/2 nematic state is characterized additionally 
by unit "pseudospin" vector cr, which changes the sign at time reversal. The additional 
mode with u^{k) has the same properties as for an easy plane Heisenberg ferromagnet, with 
the vector a playing the role of the in-plane magnetization. Thus, we arrive to the following 
non-trivial picture of the spin-3/2 nematic state: it is independently SO (3) degenerated over 
the orientation of the unit vector-director u, parallel to the long axis of the ellipsoid, and 
SO(2)-degenerated over the direction of unit pseudospin vector <?, perpendicular to u. The 
presence of this complicated spontaneous symmetry breaking leads to the appearance of 
three aforementioned Goldstone modes for this phase. 

This features of uniform nematic state gives rise for even more non-trivial properties of 
multi-sublattice state with (S) = 0. For the spin 5=1, either a two-sublattice "orthogonal 



nematic" with {u\ ■ U2) = 14j, or a "threemerized" states [17j], have been discussed. As 



we have shown, for the spin-3/2 system with bipartite lattice, two-sublattice "antinematic" 
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state, with parallel u% = U2, but antiparallel fi\ = —B2, is stable in the region 



A 2 < 0, Ai > 0. (14) 

For the rest of the parameter region, A2 < 0, Ai < 0, we have shown the presence of the 
two-sublattice antiferromagnetic state with (Si) = —(S2), |(<Si,2)| = 3/2. 

The above analysis has been done in the mean-field approximation only. A rich variety of 
"spin-liquid" states with properties governed by non-small quantum fluctuations, are known 
for one- dimensional (Id) spin-one systems (spin chains), where the mean-field approximation 
is not valid. But one can expect for the spin-3/2 system the role of such fluctuations should 
be different from those for a spin-1. In line with Haldane conjecture, for a spin 5 = 3/2 
chain one can expect a transformation of the antiferromagnetic state to a "critical" state 
with gapless elementary excitations and antiferromagnetic correlations decaying with power 
law (compare with gapped Haldane state jJLl] for integer spins 5 = 1, 2 cases). The gapless 
Luttinger liquid state for one-dimensional spin 5 = 3/2 Fermi gas has been found by Wu 12]. 
For a spin 5 = 1 chain, nematic order is probably destroyed by non-perturbative quantum 
fluctuations [5]. As has been shown above for the spin 5 = 3/2 nematic and antinematic 
states the ordering of vector-director u is accompanied by "pseudospin" a ordering. It is 
clear that the stability conditions for such states should differ from that for 5 = 1, a time- 
reversal "pure nematic" state. One can expect quite non-trivial behavior of the nematic and 
antinematic states for Id systems, but a detail investigation of these features is going far 
from the scope of this Letter. 

Thus, an isotropic magnet with spin 5 = 3/2 at L^O shows (at least, in mean-field ap- 
proximation) the presence of spin nematic and spin antinematic states with unique dynamic 
and static properties. The adequate description of an isotropic magnet with 5 = 3/2 needs 
consideration of all three possible spin invariants in the exchange Hamiltonian (TjQ). Unfor- 
tunately, we can not provide a concrete example of a crystalline magnet with spin 5 = 3/2 
and high enough biquadratic exchange. Moreover, to our knowledge, the values of bicubic 
exchange constant L, important for the presence of nematic states, has never been discussed 
for any real compounds. But the technique of ultracold gases of neutral atoms loaded into 
optical lattices gives the possibility of unprecedented control over the model parameters (ex- 
change integrals, in our consideration). Alkali atoms 132 Cs, as well as alkaline-earth atoms 
9 Be, 135 Ba, and 137 Ba, having hyperfine 3/2 spin, could be used for the study of quantum 



S 



phase transitions described above. 

To conclude, as for well-studied spin S = 1 models, spin S = 3/2 system within the 
mean-field approximation have two magnetically ordered phases, the ferromagnetic state 
and the antiferromagnetic state with maximal possible magnitude of spin on a site. As well, 
there are two different nematic states, in which the average spin equals zero. The domains 
of the stability for the nematic states separate the stability regions of the ferromagnetic and 
the antiferromagnetic phases, see Fig. 1. On the phase transition lines, Ai = or A2 = the 
symmetry is higher than the SO(3)~SU(2) rotational symmetry of the Hamiltonian. The 
qualitative differences between the nematic state for the S = 3/2 and S = 1 magnets are 
as follows. The spin S = 3/2 nematic phases display specific symmetry break with respect 
to the time reversal. It is exhibited by the properties of the quantum averages cubic over 
the spin components, which can be organized in "pseudospin" vector a. By virtue of this 
fact the spectrum of spin oscillations in the spin nematic state includes the third Goldstone 
(quasiferromagnetic) mode. The antinematic phase has the antiparallel orientations of a at 
different sublattices and possesses the following element of symmetry: the time reflection 
combined with the spatial translation on the nearest neighbor vector 5. 
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